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Summary 


A near-wall four-equation turbulence model is developed for the calculations of high-speed 
compicssible turbulent boundary layers. The four equations arc the k-e equations and the 0^-£0 
equations. These equations are used to define the turbulent diffusivities for momentum and heat 
fluxes, thus allowing the assumption of dynamical similarity between momentum and heat 
lianspoit to be lelaxed. The Favrc-averaged equations of motions are solved in conjunction with 
the four transport equations tor k, e, 0 2 and £ e . Calculations are compared with measurements 
and with another model predictions where the assumption of a constant turbulent Prandtl number is 
invoked. Compressible flat plate turbulent boundary layers with both adiabatic and constant 
temperature wall boundary conditions are considered. Cases where the free-stream Mach number 
as high as 10 and where the wall temperature ratio as low as 0.2 are calculated. Results for the 
range of low Mach numbers and temperature ratios investigated are essentially the same as those 
obtained using an identical near-wall k-£ model and the assumption of Pr^ — 0.9. One reason could 
be that the model constants used in the 0 2 and Eg equations have not been optimized to give the 
best lesults foi incompressible and compressible flows. In general, the numerical predictions are 
in very good agreement with measurements and there are significant improvements in the 
predictions of mean flow properties at high Mach numbers. Present results further show that the 
calculated Pr t for all cases investigated varies rapidly from about 0.5 at the wall to a maximum of 
approximately 1.6 in the near-wall region; however, it quickly settles to a constant value of 0.9 
beyond y w > 200. Theiefoie, the calculations lend credence to the Pr^ = 0.9 assumption invoked 
by other researchers. 
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Cr 
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model constant taken to be 0.1 

model constant taken to be 1.8 
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2 
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specific heat at constant pressure 
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model constant taken to be 0. 1 1 

model constant taken to be 0. 1 1 

model constant taken to be 0.1 1 

near-wall damping function for e equation 

near-wall damping function for Eg equation 

near-wall damping function for turbulent momentum diffusivity 

near-wall damping function for turbulent heat diffusivity 

instantaneous total enthalpy, = C p T + ^U k U k 

turbulent kinetic energy 
2 

normalized k, = k/u^ 

Mach number 

Mach number based on friction velocity, = u t /(yRT w ) 1/2 
instantaneous pressure 



p’ 

Pr 

Pr t 

R 

R t 

R e 

T 

Ui 
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I'i 

u, v 

u + 

u T 

uv + 

vB+ 

x. y 
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+ 

>'w 


Reynolds fluctuating pressure 

molecular Prandtl number 

turbulent Prandtl number 

universal gas constant 

turbulent Reynolds number, = k 2 /ve 

Reynolds number based on momentum thickness 

instantaneous temperature 

ith component of the instantaneous velocity 

instantaneous velocity components along x and y, respectively 

ith component of the Favre fluctuating velocity 

Favre fluctuating velocity components along x and y, respectively 

normalized mean U velocity, = (U)/u x 

friction velocity, = (i t w /fT W )L 2 

2 

normalized turbulent shear stress, = <uv>/u T 
normalized turbulent heat flux, = <v 0 >/U oo 0 oo 
coordinates along stream and normal directions 
normalized y coordinate, = yu x /v’ 
normalized y coordinate, = yu/v^, 
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Greek Letters 


a 

«, 

Y 

5 

£ 
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£ 

C 0 

E* 

* 

co 

E + 

C0 

e 

e 2 


thermal conductivity 
turbulent heat diffusivity 
specific heat ratio 
boundary layer thickness 
solenoidal dissipation rate of k 
dissipation rate of temperature variance 
dissipation rate defined as £ - 2"v ( 3 Vk/ 9 y) 2 
dissipation rate defined as £q - of0V02/dy)2 
dissipation rate defined as E - 2~vk/y 2 
dissipation rate defined as £y - oT02/y2 

4 

normalized dissipation rate, = £ v w /u x 
normalized dissipation rate, = £9 v^ix,.©^ 
Favre fluctuating temperature 
friction temperature 
temperature variance 
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0+ 2 

normalized temperature variance, = <0 2 >/0 oo 

0 

mean component of temperature 

K 
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fluid viscosity 

Bl 

turbulent viscosity 

V 

fluid kinematic viscosity 

V l 

turbulent kinematic viscosity, = |i t /p 

$ 

near-wall correction to £ equation 

^£0 

near-wall correction to £ e equation 

P 

instantaneous fluid density 

P’ 

Reynolds fluctuating density 

Ok 

model constant taken to be 0.75 

°e 

model constant taken to be 1.45 

O 0 2 

model constant taken to be 0.75 

a £0 

model constant taken to be 1.45 

X 

shear stress 

G>i 

fluctuating vorticity 

, Subscripts 


aw 

adiabatic wall 

r 

reference condition 

: W 

wall 

oo 

free-stream condition 

Overbars 


" 

time-averaged quantities 

Brackets 


( > 
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Introduction 


1 . 

In non-isothermal turbulent flow calculations, turbulent momentum and heat fluxes need 
modeling if the governing equations are to be closed. If, in addition, the flow is compressible, the 
modeling of these fluxes are complicated by the presence of a variable mean and fluctuating density 
in the governing equations. Conventional approach is to neglect the effects ot the fluctuating 
density and to propose models for the momentum fluxes while an additional assumption is made to 
relate the heat fluxes to the modeled momentum fluxes. Proposals for the incompressible 
momentum fluxes range from one-equation to second-order closure models (Speziale 1991). Most 
closure schemes for compressible flows invoke Morkovin's (1962) hypothesis of dynamical field 
similarity between compressible and incompressible flows. This assumption, therefore, allows the 
direct extension of incompressible models to account for compressibility effects. In addition, the 
assumption of dynamical similarity between turbulent heat and momentum transport is invoked and 
this permits the specification of a constant turbulent Prandtl number in the closure schemes (van 
Driest 1951; Anderson and Lewis 1971; Bradshaw 1977; Wilcox 1988; Speziale and Sarkar 1991; 
Aupoix and Cousteix 1991). Under these assumptions compressibility effects are accounted lor by 
the mean density alone. As a result, the ability of conventional models to reliably predict 
compressible turbulent boundary-layer flows for Mach numbers M^ > 5 has been called into 
question (Bradshaw et al. 1991; Huang et al. 1992). 

Attempts to relax some of these assumptions have been made recently. For example, 
Zhang et al. (1992) propose a compressible near-wall k-e model where all additional dilalational 
terms are systematically derived and accounted for in the governing equations. Therefore, they are 
able to assess the validity and extent of Morkovin's hypothesis. Their analysis reveals that, if the 
near-wall model is internally consistent and asymptotically correct, Morkovin's hypothesis is 
essentially valid for adiabatic wall with M m as high as 10. Consequently, the predictions in this 
Mach number range are in very good agreement with measurements. On the other hand, the model 
predictions of Cf in the case of cooled-wall boundary layers are not as good. The reason may not 



be the breakdown of Morkovin's hypothesis but rather the consequence of the assumption of a 
constant turbulent Prandtl number. The present study makes a first attempt to assess this postulate, 
and proposes to relax the assumption of a constant turbulent Prandtl number in the modeling of 
compressible turbulent boundary-layer flows. A near-wall four-equation model is suggested as an 
alternative; two equations each to model the turbulent heat and momentum fluxes. The following 
approach is adopted in the formulation of this near-wall four-equation model. 

In Zhang et al.'s (1992) near-wall k-e model, the compressible dissipation function is split 
into a solenoidal part, which is not sensitive to changes of compressibilify indicators, and a 
dilatational part, which is directly affected by these changes (Sarkar et al. 1989). This procedure, 
therefore, isolates terms in the k equation with explicit dependence on compressibility so that they 
can be modeled accordingly. An equation that governs the transport of the solenoidal dissipation 
rate with additional terms that are explicitly dependent on compressibility effects is derived 
similarly. A model (Sarkar et al. 1989) with an explicit dependence on the turbulent Mach number 
is adopted for the dilatational dissipation rate. Thus formulated, all near-wall incompressible flow 
models could be expressed in terms of the solenoidal dissipation rate and straight-forwardly 
extended to compressible flows. As a result, the incompressible equations are recovered correctly 
in the limit of constant density and vanishing turbulence Mach number. A number of near-wall 
two-equation models are available (Myong and Kasagi 1990; Deng and Piquet 1991; Karlsson et 
al. 1991; Michelassi et al. 1991; So et al. 1991a; Yang and Shih 1991). However, none is as 
widely tested for asymptotic consistency as the model of So et al. (1991a) who have validated their 
model against such benchmark data as the direct numerical simulations of channel flows (Kim et al. 
1987; Mansour et al. 1988), of flat plate boundary-layer flows (Spalart 1988) and of Couette flows 
(Lee and Kim 1991) as well as experimental measurements (Klebanoff 1955; Nishino and Kasagi 
1989). The results are in excellent agreement with data and have been reported by So et al. 
(1991b) and by Zhang and So (1991). In view of this, Zhang et al. (1992) adopt the near-wall 
model of So et al. (1991a) and extend it directly to compressible flows. Their results show that 
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compressible flat plate turbulent boundary layers can be predicted correctly up to = 10 for both 
adiabatic and cooled wall boundary conditions. Therefore, this suggests that the near-wall k-e 
model of Zhang et al. (1992) should be adopted for the momentum fluxes in the present study. 

If a constant turbulent Prandtl number is not assumed, consistent with the momentum flux 
model, a near-wall heat flux model has to be proposed. Near-wall modeling of heat fluxes is not 
as well developed; nevertheless, a second-order closure has been proposed by Lai and So (1990a) 
and a 0 2 -£ 0 model has been put forward by Nagano and Kim (1988) for non-isothermal 
incompressible flows. However, none has been formulated for compressible flows. The work of 
Zhang et al. (1992) points to the importance of having an internally consistent and asymptotically 
correct near-wall model for compressible flows. Therefore, if an incompressible near-wall model 
is to be extended to compressible boundary layers, its asymptotic behavior near a wall has to be 
analysed first. This analysis has been carried out for the 0 2 -e 0 model by Sommer et al. (1992) and 
their results show that the 0 2 -e 0 model of Nagano and Kim (1988) fails to correctly reproduce the 
asymptotic behavior of the temperature variance and its dissipation rate. Modifications along the 
line of So et al.'s (1991a) analysis of the £ equation have been proposed and a correction function 
for the £ e equation has been derived by extending the coincidence condition of Shima (1988) to the 
analysis of the £ e equation. Thus derived, the new near-wall 0 2 -£ 0 model for heat fluxes is found 
to correlate well with direct simulation data (Kim and Moin 1989; Kasagi et al. 1991) and 
experimental measurements (Johnk and Hanratly 1962; Hishida et al. 1986). In particular, the 
asymptotic near- wall behavior of the direct simulation data is reproduced correctly for both 
constant temperature and constant heat flux wall boundary conditions. This near- wall 0 2 -£ 0 model 
is extended to compressible flows in the present study by following the approach used by Zhang et 
al. (1992) to develop the compressible k-£ model. 

The four equations thus formulated are used to calculate compressible flat plate turbulent 
boundary layers with adiabatic and constant temperature wall boundary conditions. Comparisons 
with well documented experimental measurements (Femholz and Finley 1977) as well as with the 
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model calculations of Zhang et al. (1992) are carried out. Furthermore, the model is used to 
calculate the variation of skin friction coefficient with wall temperature for a cooled wall at a fixed 
free-stream Mach number and the result is compared with the van Driest II formula (Kline et al. 
1981). The validity and extent of the constant turbulent Prandtl number assumption is then 
assessed by comparing the results of the four-equation model with those of the k-e model (Zhang 
et al. 1992). 

In the following, the governing equations for compressible boundary layers are given in 
Section 2. The modeled equations for k and e as derived by Zhang et al. (1992) are presented in 
Section 3. In Section 4, the near-wall modeled equations for 0 2 and e 0 as derived by Sommer et 
al. (1992) are given together with their extension to compressible flows. Model validations and 
comparisons with measurements are discussed in Section 5. Finally, the conclusions of this study 
are presented in Section 6. 
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2 . 


The Compressible Boundary-Layer Equations 


The mean equations of motions for compressible turbulent boundary layers can be derived 
Irom the instantaneous Navier-Stokes equations by applying Favrc-averaging and then invoking 
the Prandtl boundary-layer approximations to simplify the resulting averaged equations. Favre 
decomposition is invoked for all variables except p and p where conventional Reynolds 
decomposition is assumed. In other words. 


U i = (U i ) + u i , 

(1) 

T = (0) + 0 , 

(2) 

P = ~P + p' » 

(3) 

P = P~ + p' • 

(4) 


When these decompositions are substituted into the Navier-Stokes equations and time averaging is 
applied, a set ot turbulent mean flow equations are obtained. The boundary-layer approximations 
and the assumption of negligible fluctuations in fluid properties, such as p, C p , etc, are used to 
further simplify these equations. Since the pressure field is constant for flat plate boundary layers, 
the resulting compressible turbulent boundary-layer equations can be written as (Wilcox 1988): 

Jr(p< u >) + |r(p<v>) = o . ( 5) 


P#4M SW 


dx 


dy 


dy 


dy J 


(6) 


P(U) 


d<H> , 

-dT^ pw -dT = 

d[[p Pt|d(H) 
pPr PrJ dy 




( 7 ) 
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where Cartesian x-y coordinates have been used and, consistent with conventional wisdom, the 
temperature equation is converted into the total enthalpy equation. The mean equation of state is 
assumed to be given by ’p’ = pTR <0> and Sutherland's law is used to evaluate the mean fluid 
viscosity (Zhang et al. 1992). Therefore, once ]r i and Pr t are known, (5) - (7) can be solved to 
give the velocity and temperature fields inside the boundary layers. 

In writing down these equations, gradient transport has been assumed for both the 
turbulent momentum and heat fluxes. Therefore, ifJT t is taken to be given by pTv t ), then the 
turbulent fluxes can be written as: 


- p{uv) = p v t 


m 

By 


( 8 ) 


- p (v0) = p . • (9) 

It should be pointed out that even though the equations are written in terms of a turbulent Prandtl 
number Pr t , they do not imply constant Pr t . The equations are simply written in this form for 
convenience and to comply with conventional format (Wilcox 1988; Zhang et al. 1992). Here, Pr t 
= V,/o’ l and the turbulent diffusivities are defined as: 

V t = C^k 2 /e , (10) 

a t =C x I\k[k02/ee e ]i/2 , (11) 

where the damping functions are to be defined when closure models for the momentum and heat 
fluxes arc discussed in the next two sections. In this form, Pr t varies according to the variations of 
V t and ofj. In the following, near- wall turbulence models are proposed for v t and a t so that their 
variations across the boundary layers can be determined together with other properties of the 
boundary-layer flows. 
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It should be pointed out that a rigorous derivation of (7) will give an additional turbulent 
kinetic energy term on the right hand side of (7) as pointed out by Zhang et al. (1992). In the 
present formulation, this term is omitted consistent with the work of other researcheis (Wilcox 
1988). Besides, Zhang et al. (1992) have demonstrated that the neglect of this term has a positive 
effect on the prediction of skin friction for compressible boundary-layer flows with a highly 
cooled-wall. 

The boundary conditions for <U> and <V> are no slip at the wall and <U> appioaches U M 
in the free stream. As for <H>, its free-stream value is given by = CpO**, + UJ2, while its 
wall value is taken to be either that of an adiabatic wall or a specified constant <H>. 
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3. 


Near-Wall k-£ Turbulence Model for v t 


In adopting and extending the near-wall k-e model of So et al. (1991a) to compressible 
flows, Zhang et al. (1992) made two assumptions to simplify the formulation. The first 
assumption is to split the compressible dissipation function into a solenoidal part according to the 
suggestion of Sarkar et al. (1989), so that the solenoidal dissipation is insensitive to 
compressibility effects and, therefore, approaches its incompressible limit correctly. A second 
assumption is the simultaneous neglect of fluctuating density and temperature at the wall. 
Physically, this assumption is not quite valid. However, it does permit an asymptotic analysis of 
the near-wall behavior of the compressible k and £ equations to be earned out in a manner similar 
to that proposed by Lai and So (1990b) and So et al. (1991a) for incompressible flows. 
Consequently, it is found that the near-wall k-£ model of So et al. (1991a) can be extended to 
compressible flows and the dilatational effects in the near-wall region can be accounted for by the 
varying mean density alone. Since Zhang et al. (1992) have found that the additional dilatational 
terms have little effect on compressible boundary-layer calculations in the Mach number range of 0 
< < 10 and wall temperature ratio range of 0.2 < T w /T aw < 1.0, as a first attempt, it is prudent 

to calculate turbulent heat fluxes using a near-wall model where the additional dilatational terms are 
neglected. In view of this, the near-wall k-£ model of Zhang et al. (1992) without the dilatational 
terms are adopted for the present study. The modeled k-e equations for boundary-layer flows can 
be written as: 


P(U) 


9k 

9x 


9k 

9y 


9 

(— lb] 

9k' 


(3M) 

dy . 

P + — 
\ cqj 

9y. 

+ 

9y / 


- pe 


( 12 ) 


?< u >! + p< v > 


- Pt 9e 


+ C £l |i t 


p(U) 

1 9y j 


■ c «p(f) 


+ % 


(13) 
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where £,, the near- wall function proposed for the e equation, and e, the solenoidal dissipation of k, 
are defined by 


^ = f w,2P 




(14a) 


pe = pcOiCOi 


(14b) 


Here, f w 2 is a damping function that asymptotes to one at the wall and zero far away from the 
wall. It is defined as f w 2 = e'^ Rl/64 ^“. The boundary conditions for k and £ are specified to be zero 
in the free stream. At the wall, k = 0 is assumed, while E is taken to be given by 2v w (3Vk/9y)^. 

Once k and E are known, they can be used to evaluate V t according to (10) and hence |i t = 
p~(V t ). The damping function f„ associated with v“ t is given by So et al. (1991a) as 

fp. = (1 + 3.45//RT)tanh(y + /115) . ( 15 ) 

This damping function behaves correctly as the wall is approached, i.e. f^ goes like y * 1 as y 
approaches zero. In other words, the modeled turbulent shear stress again behaves like y 3 near a 
wall similar to its exact behavior. 
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4. 


Near-Wall 0 2 -e 0 Turbulence Model for oc t 


A detailed derivation of the non-isothermal incompressible near-wall 0 2 -e 0 model has been 
given by Sommer et al. (1992). Consequently, there is no need to repeat the derivation here. 
However, some major differences between the k-e equations and the 0 2 -£ 0 equations should be 
pointed out. The first is in the modeling of the generation and destruction terms in the E 0 equation. 
Since thermal and velocity time scales are of equal importance in turbulent heat transfer, both time 
scales are used in the modeling of the generation and destruction terms. A second difference is in 
the wall boundary conditions. Constant heat flux as well as constant temperature wall boundary 
condition can be specified for non-isothermal flows. Therefore, these differences have to be taken 
into account in the derivation of the near- wall correction for the £0 equation. Sommer et al. (1992) 
incorporate the coincidence condition of Shima (1988) to treat the £q equation and derive a near- 
wall correction in a manner similar to that used by So et al. (1991a) in their derivation of %. Thus 
formulated, the 0 and £ 0 equations behave correctly as a wall is approached; at least to the lowest 
order of y. 

The near-wall 0 and £ 0 equations of Sommer et al. (1992) for a non-isothermal 
incompressible turbulent boundary-layer flow can now be written as: 


U a -^ + v de2 


3x 


3y d Y 
+ 2a t 


30' 


ikae; 

l 3y / 3y \g 0 2 3y 


‘\2 


30 

— + 2a t 

3x / 



- 2e, 


e * 



3x 


3 y 9y l 3y ) 3y (c e0 3y } 

+ C ( |2r L Pe + Qi3 £ 0 ^- - C c j4 £ 0 - 


C d i == P 0 
0 2 

Cd5^ Efl + ko > 


(16) 


(17) 
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where P is the turbulent production of k defined as P = -uv(5U /5y) and P e is the production of 
temperature variance due to mean temperature gradient and is given by P e = -[u0(5075x) + v0(5(? 
)/5y)]. The near-wall correction, ^ e0 , and the viscous dissipation of 0 , £ 0 , are defined as: 


^£0 — fw ,£0 


~ „ *2 \ 
(C d4 - 4) H e e + C d5 £ Ee - %= + (2 - Cdi - PiC|2)=p;, , 

e 2 k e 2 e 2 / 


e e = a 


50 50 
5x k 5x k 


(18a) 


(18b) 


* 

Here, P 0 = -u0 (5©/5x), where 50/5x is constant for constant wall heat flux boundary condition. 
The presence of this term in £ e0 is required in order to balance the term involving 50/5x in the 0 2 
and £ 0 equations in the near-wall region for constant wall heat flux boundary condition. It should 
be pointed out that for constant wall temperature, the term is identically zero because 50 /5x = 0. 
The damping function f w e0 is introduced to guarantee that ^ £0 vanishes far away from the wall. 
This way, the high-Reynolds-number form of the equations is recovered correctly. 

Consistent with So et al.'s (1991a) approach for the k-E model, Sommer et al. (1992) 
suggest the following form for the damping function, or f w e0 = exp[-(Re t /80) 2 ]. Once 0 2 and e 0 
are known, oc t can be determined by assuming gradient transport similar to that invoked for v t . 
However, it should be pointed out that both thermal and velocity time scales are involved in the 
transport of heat. Therefore, a t should be defined as 

a t = Cxfj t k[k9 2 /ee 0 ] 1/2 , (19) 

where the damping function f^ is given by 


fi = 


f c a 


Vr 


e t 


+ [l - exp(-y + /A + )] 2 


( 20 ) 
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The boundary conditions for 0 2 and £q are again specified as zero in the free stream and vanishing 
0 2 at the wall. As for Eq, its value at the wall is given by a v ,(d'ffi/dy) 2 . 

The incompressible near-wall 0 2 -Eq model for a t can be extended to compressible flows in 
the following manner. Again, fluctuating temperature and density are assumed to go to zero 
simultaneously at the wall and fluctuating fluid properties are neglected. Therefore, all fluid 
properties, such as p, C p , etc., can be replaced by their time-averaged values and the following 
near-wall expansions can be assumed for the fluctuating quantities. These are: 


u = a^y + a2y 2 + , 

(21a) 

v = bjy + b2y 2 + , 

(21b) 

0 = ciy + c 2 y 2 + , 

(21c) 

P' — diy + d 2 y 2 + , 

(2 Id) 


where the a's, b's, c's and d's are random functions of x, z and t. As pointed out by Bradshaw 
(1974), 0 and p' cannot go to zero simultaneously at the wall; otherwise it would lead to a zero 
wall p', which is not physically possible. In general, 0 is taken to be zero at the wall but p' is not. 
Here, p' is also assumed to be zero at the wall, however, its value away from the wall is finite. 
Therefore, this assumption represents an improvment over Morkovin’s hypothesis (1962) which 
neglects the influence of fluctuating density altogether in the whole flow. Under this assumption 
and with the help of the continuity equation for p' and (21), it can be easily shown that bj is 
identically zero irrespective of the wall thermal boundary conditions. This means that the near-wall 
asymptotic analysis of So et al. (1991a) can be used to examine the exact and modeled 
compressible 0 2 and Eq equations in the near-wall region. The result is similar to that given by 
Sommer et al. (1992) for the incompressible case. Therefore, the incompressible form of the near- 
wall 0 2 and £0 equations can be directly extended to compressible flows just like the case made for 
the k and e equations. 
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In view of this, the compressible 0 2 and eg equations can be written as: 



( 22 ) 


-/r,^ee . -* A dee 3 (--de. 


+ C dl =rp a t | 

e 2 


3 1 

f P «t 3e 0 | 

3y 1 



p a 


2 I 

+ 0(127^ p a t 

3y I k 




+ C dl ^pa t 


3y 

e a 


e 

+ c d3 ^il 


\ dx 

'w? 


V 3y 


Oi4 — P E e ' Qisr - P e e + ^e0 . 

e 2 k 


(23) 


where the near-wall correction function, % E0 , and the viscous dissipation of <0 2 >, £ 0 , are given by 


/ 


^£0 — f\V,£9 P 


*2 


(C d4 - 4) =r e e + C d5 £ e - +(2-C dl -C d2 Pr)l|P;| , 


k ' 0 2 


(24a) 


-/de 30 

e e = a( 

\ax k 3x k y 


(24b) 


r\ 

The boundary conditions for <0 > and £ 0 are the same as those specified in the incompressible 
case. However, at the wall, £ 0 is given by 0^(3^ 0 2 /3y) 2 . 
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5. 


Results and Discussion 


The governing equations (5) - (7), (12), (13), (22) and (23) are solved using the boundary- 
layer code of Anderson and Lewis (1971) with appropriate modifications made to the computer 
program. Exact boundary conditions at the wall for the turbulence quantities are used because, 
with near-wall corrections proposed in (14) and (24), the equations can be integrated directly to the 
wall and no approximations need be applied to the numerical solutions of the equations in the near- 
wall region. Three different cases are selected from the data file compiled by Fernholz and Finley 
(1977). Two cases have adiabatic wall boundary conditions and they are labeled as cases 
5301 1302 and 73050504 by Fernholz and Finley (1977). The free-stream Mach numbers of these 
two cases are 4.544 and 10.31, respectively, and the corresponding R e are 5,532 and 15,074. 
The wall temperature ratio for these two cases is O w /O r = 1.0, where 0 r is assumed to be the 
recovery temperature for adiabatic wall boundary condition and is taken to be 0 aw for cooled wall 
boundary condition. Since the fluid medium of Case 53011302 is air, Sutherland's law can be 
used to evaluate viscosity and Pr = 0.74 is specified. On the other hand, helium is used in the 
experiments of Case 73050504, therefore, Pr = 0.7. A power law as suggested by Fernholz and 
Finley (1977) is used to calculate viscosity in this case. The third case is specified by 0 w /0 aw = 
0.92, = 5.29 and R 0 = 3,939. Air is also the working fluid in this case, therefore, 

Sutherland's law can again be used to calculate viscosity and Pr = 0.74. The cases chosen span a 
wide range of and R e . In addition to comparing with measurements, the present results are 
also compared with the calculations of a near-wall k-e model where a constant Pr t is assumed. The 
model adopted is identical to solving (12) and (13) with Pr t = 0.9. Finally, another set of 
calculations is earned out so that the results can be compared with the van Driest II formula (Kline 
ct al. 1981) used to estimate the variation of Cf with 0 w /O aw for a fixed M w . This set of 
calculations is carried out at = 5.0 and Pr = 0.74. It should be pointed out that all calculations 
are earned out to the same R 0 as the experiments so that proper comparisons can be made with the 
measurements. 
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Since turbulence measurements from these three cases arc not available for comparisons, an 
alternative check on the correctness of the model is to compare the near-wall asymptotics predicted 
by the two- and four-equation models. This approach is justified because So et al. (1991a) and 
Zhang et al. (1992) have demonstrated that the near- wall asymptotics deduced from the k-e model 
are in excellent agreement with direct numerical simulations. Furthermore, Sommer et al. (1992) 
have also validated the near-wall behavior of the 0 2 and Eg equations using direct numerical 
simulation results. Therefore, a comparison of the present predictions with those of the k-e model 
can help establish the validity of the four-equation model for near-wall calculations of compressible 
flows. Once the near-wall behavior is properly established, the present calculations are used to 
assess the assumption of Pr t = 0.9 adopted in Zhang et al.'s (1992) and other researchers (Aupiox 
and Cousteix 1991; Wilcox 1988) calculations. 

According to So et al. (1991a) and Sommer et al. (1992), near a wall, the quantities k, 
<uv>, e, <0 2 >, <v0> and Eg can be expanded in terms of y. After proper normalization using 
wall variables, the expansions can be written as (So et al. 1991a; Sommer et al. 1992): 


k + = a k (yw) 2 + b k(yw) 3 + - ^ 

uv+ = a uv (yw) 3 + b uv(yw) 4 + ’ ^ 

E+ = 2a k + 4b k y+ + < 27 ) 

0 +2 = a 0 (y*) 2 + b e (yw) 3 + .... ( 28 ) 

v0 + = a v e(yw) 3 + b ve(yw) 4 + ^ 

4 = a e0 + b eeyw + •••■ (30) 


The accuracy in which the near-wall asymptotics, such as a k , a uv , a 0 , a v0 and a e0 , can be predicted 
is a measure of the correctness of the four-equation model. Furthermore, k + /e + (y^) 2 and 
0 +2 /£ 0 (y^) 2 approach exactly 0.5 and Pr, respectively, at the wall. So et al. (1991a) and Zhang et 
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al. (1992) have shown that the k-e model gives reasonable values for a k and a uv compared to direct 
numerical simulations and an exact value of 0.5 is calculated for k + /e + (y+ ) 2 at as high as 10. 
This means that the k and e equations are asymptotically correct and internally consistent. Since 
the same validation has been carried out by Sommer et al. (1992) for the 9 2 and £ e equations for 
incompressible flows, the present objective is to demonstrate that the values calculated for a 0 , a v0 
and a e0 are reasonable and that 0 +2 /£ 0 (y+) 2 is evaluated to be identical to the Pr assumed for the 
compressible calculations. Therefore, (22) and (23) can be shown to be asymptotically correct and 
internally consistent for incompressible as well as compressible flows. 

The mean velocity and temperature results are presented in Figs. 1-6 with the adiabatic 
wall cases shown in Figs. 1, 2, 4 and 5 and the cooled wall case given in Figs. 3 and 6. Three 
different ways of plotting the velocity results are presented; the conventional semi-log plot (Figs. 1 
and 3a), the linear plot (Figs. 2 and 3b) and the semi-log plot in van Driest (1951) coordinates 
(Fig. 4). On the other hand, only linear plots of the mean temperature are shown in Figs. 5 and 6 
because a friction temperature cannot be suitably defined for adiabatic wall boundary condition, 
therefore, a semi-log plot of the temperature profile is not appropriate. The rationale for presenting 
the mean velocity in these three different forms can be explained as follows. Firtsly, the 
conventional semi-log plot for compressible flows has a density effect included in the definitions of 
u + and y+ , therefore, the true velocity profile prediction cannot be directly compared with 
measurements. Secondly, errors in the predictions of the mean temperature and hence the mean 
density can occur in such a way that they tend to mask the discrepancy in the semi-log plots of the 
mean velocity. In view of this, it is also necessary to compare the mean velocity and temperature 
in linear plots so that their actual agreement with measurements can be thoroughly analysed. 
Thirdly, van Driest suggests stretching u + further by a density ratio so that a new can be defined 
as 
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With this new coordinate, the compressible law-of-the-wall as deduced by van Driest (1951) and 
simplified by Bradshaw (1977) can be written as 

= oir lny " + 5,2 + 95M * ‘ (32) 

This form differs from the conventional law-of-the wall which is given by 

u+ = oir lny ” + B ’ (33) 

where the constant B is a function of Mach number for compressible flat plate turbulent boundary 
layers. When the velocity results are plotted in terms of u + and u^, the validity of (32) and (33) 
can be evaluated. Thus presented, the mean properties can be thoroughly compared and their 
agreement or lack thereof with measurements and other model calculations can be analysed. 
Finally, the calculated Cf for the different cases is tabulated in Table 1 for comparison with data. 

It can be seen from these plots that the model calculations are in good agreement with 
measurements and the predictions of the k-e model. Furthermore, these results are essentially 
identical at low Mach numbers for both adiabatic (Figs, la, 2a and 5a) and cooled wall (Figs. 3 
and 6) conditions. The only difference appears to occur in the case of M^ = 10.31, where the 
four-equation model gives a significant improvement in the predictions of both mean velocity and 
temperature (Figs, lb, 2b and 5b) compared to those given by the k-e model. This improvement is 
due to a better estimate of the turbulent Prandtl number near the wall. More will be said about this 
when the turbulence properties in the near-wall region are examined. Zhang et al. (1992) have 
demonstrated that the calculated u + can be described fairly well by the conventional law-of-the-wall 
(33) and the constant B thus deduced is approximately 4.7 for the three cases examined. The 
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present results are in agreement with their conclusion. Therefore, the variable turbulent Prandtl 
number formulation has little or no effect on the log region of the boundary-layer flow. Plots of 
the mean velocities in van Driest coordinates for the case with = 4.544 are shown in Fig. 4 
together with a plot of (32). A line parallel to (32) can be drawn through some of the data points; 
however, the intercept thus deduced is different from that given in (32). On the other hand, the 
calculated profiles from the two different models are in very good agreement with data over a much 
wider range of y^, and the slope of the log-law thus determined yields a von Karman constant k = 
0.35 which is significantly smaller than a value of 0.41 quoted by Bradshaw (1977). Finally, the 
calculated C f 's are compared with data in Table 1 and, in general, the four-equation model gives a 
slight improvement over that of the k-e model. 

The turbulence properties in the near- wall region are plotted in Figs. 7-12. Only the 
profiles of k + , e + , uv + and v0 + are compared. The results for the adiabatic wall cases are 
presented in Figs. 7, 8, 10 and 1 1 while those for the cooled wall case are shown in Figs. 9 and 
12. At low Mach numbers, the predictions of these properties by the two different models are 
essentially identical. This conclusion applies to both adiabatic (Figs. 7a, 8a, 10a and 11) and 
cooled wall (Figs. 9 and 12) boundary conditions. The only difference comes in the predictions of 
the case where = 10.31. In general, the present calculations of k + , e + and uv + are slightly 
lower than those predicted by the k-e model. One of the reason is the variable turbulent Prandtl 
number. Plots of Pr t across the inner region of the boundary layers are shown in Figs. 13 and 14. 
The turbulent Prandtl number is seen to vary rapidly in the region very near the wall. It increases 
from a wall value of about 0.5 to a maximum of approximately 1.6 and then decreases to about 
0.75 before settling back to a value of 0.9 at y^ = 200. Thereafter, Pr t remains fairly constant at 
0.9. This shows that all variations of Pr t occur in the region, 0 < y^, < 200. Consequently, it is 
not surprising to find that differences in k + , e + and uv + between the two model calculations also 
lake place in this region for the case where = 10.31. In conclusion, it can be said that variable 


18 


Pr, effect on the calculated properties is small and this observation is also supported by a 
comparison of the near-wall asymptotics which are tabulated in Table 2. 

In general, the a k values calculated by the four-equation model are slightly lower than those 
deduced from the k-E model and the differences are insignificant. The other quantities calculated 
by the two different models are essentially the same. The four-equation model again yields a 
k + /e + (yw) 2 = 0-5 for all three cases considered, thus showing that it is asymptotically correct and 
internally consistent as far as the k and £ equations are concerned even though they are coupled to 
the 0 and Eq equations through the mean velocity and mean density fields. Since a friction 
temperature cannot be defined for adiabatic wall boundary condition, the free-stream temperature 
has to be used to normalize <0 2 > and Eg. As a result, their calculated near-wall asymptotics are at 
least one order of magnitude smaller than those for k and E. The calculated a v9 from the two 
different models are of the same order with those deduced from the four-equation model generally 
higher than those determined from the k-£ model. Furthermore, the 0 and Eg equations are also 
asymptotically correct and internally consistent because the calculations of 0 +2 /e^ (y^,) 2 are identical 
to the Pr assumed for each case. Once the reliability of the near-wall asymptotics has been 
established, the calculated variations of Pr, become credible. For all three cases considered, the Pr, 
calculated is constant beyond y* = 200 and its value is 0.9. This result lends credence to the Pr, = 
0.9 assumption invoked by other researchers (Wilcox 1988; Aupoix and Cousteix 1991; Zhang et 
al. 1992) and suggests that it is essentially valid for the range of M TC and 0^0,^ considered. 

Zhang et al. (1992) have shown that the k-E model gives excellent prediction of the 
variation of Cf/(Cf), with in the range 0 < M m < 10 for adiabatic wall boundary condition 
compared to the van Driest II formula (Kline et al. 1981). Since the present model is in good 
agreement with the k-E model in its prediction of C f for adiabatic wall boundary condition, there is 
no further need to verify the validity of the four-equation model for its ability to calculate correctly 
the variation of Cf/(C f )i with in the range 0 < M m < 10. On the other hand, its ability to predict 
the variation of Cf/(Cf), with 0 w /0 aw for a fixed M w has to be verified. A comparison of the 


19 



present calculations at a constant M„ = 5.0 with the results deduced from van Driest II formula 
(Kline et al. 1981) and the k-e model is tabulated in Table 3. As suggested by Kline et al. (1981), 
the calculations are earned out to Rq = 10 4 and (Cf), = 2.70xl0' 3 is assumed. The results reveal 
that there are essentially very little difference between the present calculations and those obtained 
from the k-e model. Both sets of predictions are in good agreement with the van Driest II formula. 
Therefore, the ability of the four-equation model to calculate Cf correctly for cooled wall boundary 
coditions is established. 
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6 . 


Conclusions 


A near-wall four-equation turbulence model has been developed for the calculations of 
compressible flat plate turbulent boundary layers with constant heat flux and constant temperature 
wall boundary conditions. The four equations consist of the transport equations for k, e, 9 and 
e e . These equations are modified for near-wall flow calculations so that they can be integrated 
directly to the wall and the exact boundary conditions at the wall can be satisfied. The 

modifications of So et al. (1991a) for the k and £ equations are adopted and extended directly to 

— — ■ 

compressible flows. Similar modifications for the 0 and £ e equations have been carried out for 
the incompressible form of these equations (Sommer et al. 1992) and they are extended to 
compressible flows by invoking Morkovin's (1962) hypothesis. Thus formulated, the four- 
equation model is internally consistent and asymptotically correct near a wall, and there is no need 
to assume a constant turbulent Prandtl number because the turbulent heat flux can be estimated 
from a knowledge of 0 and £ 0 by assuming gradient heat transport. 

The near-wall four-equation turbulence model is used to calculate compressible flat plate 
turbulent boundary layers with free-stream Mach numbers as high as 10 and with adiabatic and 
cooled wall boundary conditions. The calculations are compared with measurements and with the 
predictions qf a near-wall k-£ model where Pr t = 0.9 is assumed. Three cases have been 
calculated; two with adiabatic wall boundary condition and one with constant wall temperature. 
The free-stream Mach numbers for these cases vary from a low of 4.544 to a high of 10.31 and the 
wall temperature ratio varies from 0.2 to 1. Good agreement with measurements and the k-£ model 
calculations is obtained. In general, a variable turbulent Prandtl number formulation improves the 
calculated properties; particularly for high free-stream Mach numbers. At M^ = 10.31, there arc 
significant improvements in the predictions of the mean velocity and mean temperature compared to 
the k-£ model. The turbulent Prandtl number thus calculated has a wall value of about 0.5 for all 
cases considered. It increases sharply to approximately 1.6 away from the wall and then settles 
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down to 0.9 at a y^, = 200. These results, therefore, verify that the Pr t = 0.9 assumption invoked 
by past researchers is valid for compressible flat plate turbulent boundary layers for the range of 
Moo and O w /0 aw examined in this study. 
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Table 1. Comparison of calculated and measured Cf. 


J 


Source 

CfX 10 3 

i 

1 

M« = 4.544 

M m = 10.31 

= 5.29 

! 

Data (Fernholz and Finley 
1977) 

1.26 

0.24 

1.31 

k-e model (Zhang et al. 

1.320 

0.220 

1.262 


1992) 




t 

Four-equation model 

1.35 

0.24 

1.28 

z 
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Table 2. 


Comparisons of the near- wall asymptotics deduced from the two different models. 


Near- wall 
asymptotics 

k-e model (Zhang et al. 1992) 

Four-equation model 

= 4.544 

M oc = 10.31 

- 5.29 

= 4.544 

M co = 10.31 

Moo = 5.29 

a k 

0.0836 

0.0771 

0.0788 

0.0824 

0.0739 

0.0772 

a uv x 10 4 

6.76 

6.74 

6.14 

6.75 

6.30 

6.10 

k + /e + (Yw) 2 

0.50 

0.50 

0.50 

0.50 

0.50 

0.50 

a e 

- 

- 

- 

0.00310 

0.23415 

0.00547 

a e9 

- 

- 

- 

0.00419 

0.33450 

0.00739 

a v0 x 10 5 

- 0.1 18 

- 1.310 

- 0.059 

- 0.20 

- 1.69 

- 0.10 

e +2 /4(yw) 2 

- 

- 

” 

0.74 

0.70 

0.74 
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Table 3. Comparisons of calculated C f /(Cf)j for cooled wall boundary condition at M 


®v/®aw 

Cf/(C f )i 

van Driest II 
(Kline et al.1981) 

k-e model (Zhang et 
al. 1992) 

Four-equation 

model 

0.2 

0.58 ± 0.058 

0.52 

0.52 

0.4 

0.49 ± 0.049 

0.47 

0.49 

0.6 

0.43 ± 0.043 

0.45 

0.46 

0.8 

0.39 ± 0.039 

0.42 

0.42 






1.0 


0.35 ± 0.035 


0.38 


0.38 
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o.o 0.1 



o.o 0.1 


0.5 0.6 0.7 0.8 


Figure 2. Comparison of the mean velocity in linear plots for adiabatic 
wall boundary condition: (a) M«> = 4.544, (b) M« = 10.31. 





Figure 3. Comparison of the mean velocity for cooled wall condition at 
Moo = 5.29: (a) semi-log plot, (b) linear plots. 
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Comparison of the mean temperature in linear plots for adiabatic 
wall boundary condition: (a) M» = 4.544, (b) M„o = 10.31. 







Figure 6. Comparison of the mean temperature in linear plots for cooled 
wall boundary condition at M«, = 5.29. 
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Figure 12. Comparison of (a) <uv> and (b) <v0> for cooled wall boundary 
condition at M»d = 5.29. 
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Figure 13. 


Variation of Pr t across the boundary layer for adiabatic wall 
boundary condition: (a) M» = 4.544, (b) Moo = 10.31. 







Figure 14. Variation of Pr t across the boundary layer for cooled wall 
boundary condition at M* = 5.29. 
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